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Abstract 

We investigate the thermal Casimir interaction of a suspended graphene described by the Dirac 
model with a plate made of dielectric or metallic materials. The reflection coefficients on graphene 
expressed in terms of a temperature-dependent polarization tensor are used. We demonstrate that 
for a graphene with nonzero mass gap parameter the Casimir free energy remains nearly constant 
(and the thermal correction negligibly small) over some temperature interval. For the interaction 
of graphene with metallic plate, the free energy is nearly the same, irrespective of whether the 
metal is nonmagnetic or magnetic and whether it is described using the Drude- or plasma-model 
approaches. The free energy computed using the Dirac model was compared with that computed 
using the hydrodynamic model of graphene and big differences accessible for experimental obser- 
vation have been found. For dielectric and nonmagnetic metallic plates described by the Drude 
model these differences vanish with increasing temperature (separation). However, for nonmagnetic 
metals described by the plasma model and for magnetic metals, a severe dependence on the chosen 
theoretical description of graphene remains even at high temperature. In all cases the analytic 
asymptotic expressions for the free energy at high temperature are obtained and found in a very 
good agreement with the results of numerical computations. 

PACS numbers: 78.67.Wj, 42.50.Lc, 65.80.Ck, 12.20.-m 
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I. INTRODUCTION 



It has been known that nanostructures based on carbon possess unique mechanical, elec- 
trical and optical properties.- Among them particular attention has been given to graphene, 
a two-dimensional sheet of carbon atoms arranged in a hexagonal structure with low-energy 
electronic excitations described by the Dirac equation.- At the present time suspended 
graphene membranes up to 55 /zm diameter are produced. 3 This makes possible to investigate 
the van der Waals and Casimir interaction between graphene and different material struc- 
tures, such as another sheet of graphene, atoms, molecules, dielectric and metallic plates, 
spheres etc. Respective theoretical investigations were performed using the phenomenolog- 
ical density-functional methods,- - - second-order perturbation theory^ and, for multilayered 
carbon nanostructures, using the Lifshitz theory.— 

To apply the Lifshitz theory, one needs the reflection coefficients on graphene over a 
wide frequency region. It is customary to express the reflection coefficients in terms of the 
frequency- dependent dielectric permittivity, a concept which is not well defined for one- 
atom-thick carbon nanostructures. Because of this, two models for the reflection coefficients 
on graphene with no use of dielectric permittivity were proposed, the hydrodynamic one^^- 
and the Dirac one.— ^ In the framework of the hydrodynamic model, graphene is considered 
as an infinitesimally thin positively charged sheet, carrying a homogeneous fluid with some 
mass and negative charge densities. The hydrodynamic model was applied for calculation of 
the Casimir and Casimir-Polder interactions.— ~— In the framework of the Dirac model, it is 
taken into account that for energies below a few eV the dispersion relation for quasiparticles 
in graphene is linear with respect to the momentum, whereas it is quadratic in the hydrody- 
namic model. The reflection coefficients of graphene at zero and nonzero temperature were 
found in Refs.— ancU^, respectively. Some calculations of the Casimir-Polder graphene-atom 
interaction were performed at zero temperature using both the hydrodynamic and Dirac 
models.— 1 ^ It should be remembered, however, that both these models are only approxima- 
tions. As was already mentioned, the hydrodynamic model disregards the Dirac character of 
charge carriers in graphene. As to the Dirac model, it extends the linear dispersion relation 
for quasiparticles to any energy, whereas this property applies only at low energies. Because 
of this, one can conclude that in calculations of the Casimir and Casimir-Polder forces the 
Dirac model of graphene should be applicable at large separations between the test bodies, 
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whereas the hydro dynamic model might work at short separation distances. 

There is also another approach to the application of the Lifshitz theory to graphene. 
In this approach graphene is characterized by a spatially nonlocal dielectric permittivity 
depending both on the frequency and the wave vector. Different authors express such a 
dielectric permittivity either through the polarizability of one graphene layer— or in the 
random phase approximation.— 

It is well known that thermal Casimir effect is a subject of debate and there are differ- 
ent theoretical approaches to its description.— _ — For two graphene sheets in a nonretarded 
regime it was found^ that a relatively large thermal correction to the Casimir force at room 
temperature arises at short separation distances of tens of nanometers. This conclusion 
was qualitatively confirmed,-^ using the Dirac model of graphene, with an alternative ex- 
planation. The reason for the origin of large thermal correction for graphene is that the 
contribution of all terms with nonzero Matsubara frequencies at room temperature becomes 
small in comparison with the zero-frequency term even at short separations. It was also 
found 28 that in graphene-atom Casimir-Polder interaction the thermal effect depends cru- 
cially on the magnitude of a mass gap parameter in the Dirac model. Specifically, for a 
nonzero gap there exists an interval of temperatures (separations) where the thermal cor- 
rection remains small with increasing temperature (separation). The possibility of large 
thermal correction at short separations links the Dirac model of graphene to the Drude 
model used in the literature to describe the Casimir effect between real metals.— 1 ^ Because 
of this, it is of much interest to investigate the thermal Casimir effect in the interaction of 
graphene with real material bodies made of different materials. 

In this paper, we calculate the free energy of the Casimir interaction between a sus- 
pended graphene membrane described by the Dirac model and dielectric (silicon, sapphire) 
or metallic (Au, Ni) plates. In so doing materials of the plate are described by realistic 
dielectric permittivities taking into account the interband transitions of core electrons. We 
demonstrate that, similar to graphene-atom interaction, the behavior of thermal correction 
crucially depends on the mass gap parameter A of the Dirac model. Note that although the 
Dirac-type excitations in pristine graphene are gapless, the influence of electron-electron in- 
teraction, substrates, defects of structure, and other effects leads to a nonzero mass gap.->22r— 
Specifically, we show that larger is the magnitude of mass gap parameter, wider is the sep- 
aration (temperature) interval, where the thermal correction remains small with increasing 
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separation (temperature). For a metallic plate interacting with graphene, we perform all 
calculations using both the Drude- and plasma-model approaches to the dielectric permit- 
tivity of metal. In the case when graphene described by the Dirac model interacts with a 
metallic plate (either nonmagnetic or magnetic) the calculation results obtained using both 
approaches nearly coincide and do not depend on the magnetic properties. In contrast to the 
case of two Drude metals, the thermal correction for a graphene-metal interaction has the 
same sign as the interaction energy at zero temperature, i.e., the magnitude of the Casimir 
free energy increases with increasing temperature. Note that all results obtained for a free 
energy in graphene-plate geometry can be reformulated as the Casimir force between a ma- 
terial sphere and a graphene sheet using the proximity force approximation (PFA).— As was 
recently shown,— ~— the error arising from the use of PFA is less than the ratio of separation 
distance to sphere radius. 

In this paper we also compare the predictions of the Dirac model for the thermal Casimir 
effect with respective predictions of the hydrodynamic model and discuss the application 
region of each. Specifically, it is shown that the hydrodynamic and Dirac models of graphene 
lead to different results at short separations and to nearly coinciding results at large sep- 
arations for the Casimir free energy of graphene interacting at room temperature with a 
dielectric plate or with a nonmagnetic metallic plate described by the Drude model. For 
a nonmagnetic metallic plate described by the plasma model or for a magnetic plate the 
predictions of the hydrodynamic and Dirac models of graphene are significantly different at 
all separations considered and can be discriminated experimentally. At large separations the 
asymptotic expressions for the Casimir free energy of graphene-plate interaction are derived 
and compared with the computational results. 

The paper is organized as follows. In Sec. II we present the reflection coefficients of the 
electromagnetic oscillations on graphene in the Dirac and hydrodynamic models. Section III 
is devoted to the thermal Casimir interaction of graphene described by the Dirac model with 
a dielectric plate made of silicon or sapphire. Similar results for graphene interacting with 
a metallic plate made of Au and Ni are presented in Sec. IV. In Sec. V the theoretical pre- 
dictions following from the Dirac and hydrodynamic models are compared. Our conclusions 
and discussions are contained in Sec. VI. 



4 



II. REFLECTION COEFFICIENTS ON GRAPHENE 



Here, we briefly present the Lifshitz formula for the free energy of graphene interacting 
with a material plate and respective reflection coefficients on graphene derived using the 
Dirac and hydrodynamic models. It is supposed that the suspended graphene is at a sep- 
aration a from the thick plate (semispace) at thermal equilibrium at temperature T. The 
material of a plate is described by the frequency-dependent dielectric permittivity s(ou) and 
magnetic permeability /x(u;). The Casimir free energy per unit area J 7 is given by the Lifshitz 
formula. 33 - For simplicity in computations, we express it in terms of dimensionless variables 
as follows: 



^' t )=^Se7 ydy 

1=0 J 

+ In [l - r g(iO, y)rg(tCi, v)**] } • (1) 

Here, ks is the Boltzmann constant, Q are the dimensionless Matsubara frequencies con- 
nected with the dimensional ones, £j = 2irkBTl/h, by the equality Q = £i/cu c where 

= c/(2a). The dimensionless wave vector variable y is connected with the magni- 
tude of the projection of the wave vector on the plane of a plate, k±, by the equality 
y = 2a(kj_ + Cf/ c 2 ) 1 / 2 . The reflection coefficients on graphene, ^tm(te)' anc ^ on a pl a te, 
r TM(TE)' are ^ or t wo independent polarizations of the electromagnetic field, transverse mag- 
netic (TM) and transverse electric (TE). They are taken at imaginary frequencies. The 
prime near the summation sign means that the term with I = is taken with a factor 1/2. 

The reflection coefficients on graphene in the Dirac model were expressed in terms of 
the components of the polarization tensor in three-dimensional space-time in the following 
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Here, the dimensionless components of the polarization tensor are defined as 

noQ,tr — -JTi±00,tr (3) 



and trace stands for the sum of spatial components Hi and IT|. 

At nonzero temperature the explicit expression for the polarization tensor in the Dirac 
model with arbitrary mass gap parameter A and chemical potential /i was found in Ref.—. 
Here we consider the case of undoped graphene and put \i = 0. Then in terms of our 
dimensionless variables the result of Ref.— for the 00-component of the polarization tensor 
takes the following equivalent form:— 
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In this expression, a = e 2 /(hc) in the fine-structure constant (e is the electron charge), 
A = A/ \TvjJ c ) is the dimensionless mass gap parameter, the dimensionless Fermi velocity 
is v F = vp/c ~ 1/300, and the dimensionless variable r is defined as r = 2nT/T e Q = 
AnaksT '/ '(he). Equation (jl]) also contains two dimensionless functions defined by 
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The result of Ref.— for the trace of the polarization tensor in terms of our dimensionless 
variables is given by2£ 
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The properties of the reflection coefficients ([2]) with the polarization tensor (j3J) and ([6]) were 
studied previously.— 1 ^ 

In the framework of the hydrodynamic model discussed in Sec. I the reflection coefficients 
on graphene take the more simple form^i^^^ 

Ky + Cf 
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Here, the dimensionless characteristic wave number of the graphene is defined as K = 2aK, 
where the dimensional wave number is 
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K = 2vr — - = 6.75 x 10 5 m _1 . (8) 

mc 2 

In Eq. flSJ n is the number of 7r-electrons per unit area, m is the electron mass. Note that 
the parameter K of the hydrodynamic model does not depend on temperature. Thus, the 
reflection coefficients (|7j) and the free energy ([TJ) in the case of hydrodynamic model depend 
on the temperature only through the Matsubara frequencies. A different situation arises 
in the Dirac model. Here, the reflection coefficients (121) depend on T not only through the 
Matsubara frequencies, but also through the components of the polarization tensor. This 
links the Dirac model of graphene to the Drude-model approach to the thermal Casimir 
force between real metals (because the relaxation parameter of the Drude model and, thus, 
the reflection coefficients are the explicit functions of temperature). In this respect the use 
of the hydrodynamic model of graphene in calculations of the Casimir force is analogous to 
the plasma-model approach where the reflection coefficients depend on T only through the 
Matsubara frequencies. 

As to the reflection coefficient of electromagnetic oscillations on thick metallic plate 
(semispace), they have the standard form 33 
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Here, both the dielectric permittivity E\ = e{iu c Cj) and the magnetic permeability \L\ 
n(iu) c (i) are calculated along the imaginary frequency axis. 



III. THERMAL INTERACTION OF GRAPHENE DESCRIBED BY THE DIRAC 
MODEL WITH DIELECTRIC PLATE 

Here, we investigate dependences of the free energy as functions of temperature and 
separation in the Casimir interaction of a suspended graphene sheet with the dielectric 
plate made of different dielectrics. As the most typical dielectric materials we consider 
silicon (Si) and sapphire (AI2O3) which possess relatively large values of the static dielectric 
permittivity [e(0) = 11.67 and 10.102, respectively] but quite different behaviors along the 
axis of imaginary frequencies (for Si it is caused by electronic polarization, whereas for Al 2 03 
by both electronic and ionic polarizations). 

A. Free energy as a function of temperature 

We begin with the free energy of interaction of graphene described by the Dirac model 
with Si plate. The dielectric permittivity of Si along the imaginary frequency axis was 
obtained^I^ by means of the Kramers-Kronig relations from the tabulated optical data for 
the complex index of refraction.— It is assumed to be temperature- independent. Compu- 
tations were performed by Eq. ([!]), where the reflection coefficients on graphene are given 
by Eqs. fl2]), (0J and flS]), and on silicon by Eq. (jUJ) with \i\ — 1 and Si specified above, 
over the temperature interval from to 300 K at two separation distances a = 100 nm and 
a = 1 /im. For the mass gap parameter of the Dirac model of graphene only the upper bound 
is known.— i22a22 For a suspended graphene we choose the realistic upper bound A < 0.1 eV. 
For a graphene deposited on substrate, A can be several times larger.— Taking this upper 
bound into account, we perform all computations for A = 0.1, 0.05, 0.01, and OeV. 

The computational results for the free energy per unit area are presented in Fig. [T] as 
functions of temperature, where the lines from top to bottom correspond to A decreasing 
from 0.1 eV to OeV, (a) at the separation a = 100 nm and (b) at a = 1 /im. As can be seen 
in Fig. QJa,b), for each mass gap A/0 there exists the temperature interval where the free 
energy remains nearly constant with increasing temperature. The width of these intervals 
quickly decreases with decreasing A. Thus, at a = 100 nm and A = 0.1 eV the free energy is 
nearly constant up to T = 120 K. At the same separation, but with A = 0.05 eV, the same 
property holds only up to T = 70 K. The widths of intervals, where the free energy remains 



nearly constant with increasing T, are also narrowed with the increase of separation [see 
Fig. QJb)]. Note that computations performed for A < 0.001 eV lead to nearly the same 
numerical results as for A = [the lowest lines in Fig. []Ja,b)]. 

Figure [T] suggests that within the temperature interval, where the free energy is nearly 
flat, the thermal correction to the Casimir energy at zero temperature should be relatively 
small. We confirm this conclusion by the direct computation of the thermal correction to 
the Casimir energy defined as 

A T F(a,T) = T(a,T) - T(a,0) (10) 

for the same values of parameters, as in Fig. [TJ 

The computational results for the thermal correction as a function of temperature are 
presented in Fig. [2]at separations (a) a = 100 nm and (b) a = 1 /im. The lines labeled 1, 2, 3, 
and 4 correspond to the values of mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, 
respectively. As expected, for each line with A^0 there is some interval where the thermal 
correction remains nearly zero. These intervals are just the same where the Casimir free 
energy in Fig. [1] remains nearly flat. It is interesting to note that the thermal corrections are 
monotonously decreasing functions of temperature and have the same negative sign as the 
free energy. This contrasts with the Drude model approach to the Casimir force between 
real metals where the thermal correction over a wide temperature interval is positive making 
the free energy the nonmonotonous function of temperature.— 1 ^ 

It is instructive also to compute the relative thermal correction to the Casimir energy 
defined as 

A T F(a,T) 

5 T F{a, T) = . 11 

^{a, 0) 

The computational results are presented in Fig. [3] in percent as functions of temperature at 
separations distances (a) a = 100 nm and (b) a = 1 fim. The lines 1, 2, 3, and 4 are labeled 
in the same way as in Fig. |2j As can be seen in Fig. [3j all the relative thermal corrections 
are monotonously increasing functions of temperature whose character depends crucially on 
the magnitude of a mass gap parameter A. For example, at a = 100 nm [Fig. [3]^a)] the 
relative thermal correction of line 1 (A = 0.1 eV) is nearly equal to zero below T = 120 K 
but achieves 120.44% at T = 300 K. To compare, at the same separation the relative thermal 
correction of line 4 (A < 0.001 eV), which is much larger than the corrections of lines 1-3 at 
low temperatures, achieves only 59.40% at T = 300 K. It is seen that at room temperature 



the relative thermal correction in the Casimir interaction of graphene with Si is rather 
large even at relatively short separations. From Fig. [3](b) it follows that at a = 1 /zm the 
relative thermal correction at room temperature exceeds 4000%, i.e., the absolute thermal 
correction exceeds the Casimir energy at zero temperature by a factor of 40. This makes 
graphene interesting for experimental investigation of thermal effects in the Casimir force. 

Next we consider the interaction of graphene described by the Dirac model with sapphire 
plate. The dielectric permittivity of sapphire along the imaginary frequency axis is well 
described 40 in the Ninham-Parsegian approximation and was already used^ in computations 
of the Casimir force. Computations of the free energy as a function of temperature were 
performed in the same way as for silicon. The computational results for a = 100 nm are 
presented in Fig. HJa) by the four lines from top to bottom for A = 0.1, 0.05, 0.01, and 
< 0.001 eV, respectively. Similar to Fig. [Ha), there are temperature intervals where the 
free energy remains nearly constant. For sapphire, however, the respective magnitudes 
of the free energy are smaller than for a silicon. Skipping the computational results at 
a = 1 /im [which are similar to those presented in Fig. DQ[b)], we present in Fig. H^b) the 
more detailed computational results for the free energy in the temperature interval from 
OK to 200 K, where the lines from top to bottom correspond to A = 0.1, 0.09, 0.08, 0.07, 
0.06, and 0.05 eV, respectively. Keeping in mind that the exact value of A is not known, 
computations of this kind can be useful for the determination of A from the comparison 
between experiment and theory. For sapphire, the computational results for the absolute 
and relative thermal corrections are similar to those for silicon (see Figs. [2] and E]). Because 
of this we do not present them here. 

B. Free energy as a function of separation 

Here, we present the computational results for the interaction of graphene described by 
the Dirac model with dielectric plates as a function of separation. The same equations and 
dielectric functions, as in Sec. IIIA, are used. Taking into account that the Casimir free 
energy strongly depends on separation, we normalize the results obtained on the Casimir 
energy per unit area between two parallel plates made of ideal metal 
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In Fig. EI a) the quantity J 7 / Ec at T = 77 K is plotted as a function of separation over 
the region from 50 nm to 5 fim. The lines from bottom to top correspond to the mass gap 
parameter equal to 0.1, 0.05, and < 0.01 eV, respectively. In all cases the magnitude of the 
free energy decreases monotonously with the increase of separation (the increase of T /Eq 
is Fig. is explained by the fact that \Ec\ decreases with separation faster than \J Z \). In 
Fig. Elb) similar results are shown at T = 300 K. Here, to avoid an overlap of the lines, 
we consider a more narrow separation region from 50 nm to 1 yum. The bottom and top 
lines correspond to A = 0.1 and < 0.01 eV, respectively. As can be seen from Fig. Elb), 
at T = 300 K the dependence of the Casimir free energy on the mass gap parameter of 
the Dirac model is very weak. Moreover, at separations above 50 nm the magnitude of the 
free energy T decreases with separation as a -2 , so that the quantity J- ' /Ec is nearly linear 
function of separation. 

Similar computations were performed for a graphene interacting with a sapphire plate. 
In Fig. Ela) the computational results for J- '/Ec as a function of separation at T = 77 K are 
presented, where the lines from bottom to top correspond to A = 0.1, 0.05, and < 0.01 eV, 
respectively. These results are similar to those presented in Fig. Ela) for a silicon plate. 
Again, at T = 77 K the Casimir free energy is strongly affected by the value of A. To 
illustrate this in more detail, in Fig. Elb) we plot J^/Ec over the separation region from 
50 nm to 2 fim where the lines from bottom to top correspond to A = 0.1, 0.09, 0.08, 0.07, 
0.06, and 0.05 eV. At T = 300 K, similar to Fig. Elb), the dependence of the free energy on 
A becomes very weak and one obtains J^/Ec increasing nearly linear with the increase of a. 

C. Asymptotic behavior at high temperature 

The asymptotic behavior of the Casimir free energy at high temperature (or, equivalently, 
at large separations) can be obtained by considering the zero-frequency contribution to the 
Lifshitz formula ([1]). This corresponds to large values of the dimensionless parameter r 
introduced after Eq. @. By putting I = 0, Co — in Eq. (0J one arrives at 
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In the limit of high temperature we assume that irfp/r <C 1. In this case 
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From Eq. (j2J) the TM reflection coefficient on graphene at zero Matsubara frequency is given 
by 
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Taking into account that in accordance with Eq. (jUJ) for dielectric materials r^(0, ?/) = 0, 
we obtain from Eqs. ([1]) and (ITBl the following asymptotic expression for the Casimir free 
energy at large r: 
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where the TM reflection coefficient of the dielectric plate at zero Matsubara frequency 
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In view of the fact that t/(ttvf) S> 1, one obtains 



Performing the integration in Eq. 
the Casimir free energy: 
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where Li n (z) is the polylogarithm function. 

The application region of Eq. (123]) depends on the specific values of parameters. Thus, 
for Si at T = 300 K Eq. ( 12"B"|) leads to less than 1% errors in the values of the free energy, 
as compared with the results of numerical computations, at a > 500 nm for graphene with 
A = 0.1 eV and at a > 1 /j,m for graphene with A = 0.01 eV. For Si at T = 77 K and 
graphene with A = 0.1 eV Eq. (123]) is not yet applicable at a = 5 /mi and for graphene with 
A = 0.01 eV works well for a > 4 fim. 

It is interesting to compare the asymptotic expression (123]) with other results obtained 
in the literature. Thus, using the nonlocal dielectric function in the random phase approx- 
imation, the free energy of graphene interacting with a dielectric substrate (SiC^) at large 
separations was found^ to decrease as a" 3 . This is not in accordance with the main term 
of our result ( 123]) which demonstrates the classical limit, as is expected at large a. Note 
that another work— models the dielectric properties of graphene by the Drude-type function 
and arrives at the a~ 2 scaling for graphene-graphene interaction which satisfies the classical 
limit. 



IV. THERMAL INTERACTION OF GRAPHENE DESCRIBED BY THE DIRAC 
MODEL WITH METALLIC PLATE 

The case of graphene interacting with metallic plate is of special interest. As was men- 
tioned in Sec. I, there are two different theoretical approaches to the description of Casimir 
effect between real-metal plates. The Drude model approach takes into account the relax- 
ation properties of conduction electrons. In the framework of this approach, the imaginary 
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part of the dielectric permittivity of the Drude model is used to extrapolate Ime opt (u), 
obtained from the measured optical data, to zero frequency. By contrast, the plasma model 
approach disregards relaxation processes and extrapolates e opt (u) to zero frequency by means 
of the simple plasma model. Although the Drude model approach may seem preferable, as it 
takes into account some really existing property of metals, the experimental situation more 
likely favors the plasma model approach. In a series of precise independent measurements 
performed by the two experimental groups^r— the Drude model approach was excluded at 
a high confidence level (several experiments^ - — also excluded the influence of free charge 
carriers, that are present in dielectric materials at room temperature, on the Casimir force). 
The two experiments that support the Drude model approach 5 ^ 1 ^ are not independent mea- 
surements of the Casimir force; they are based on fitting procedures between measured data 
for the total force and theoretical predictions using hypothetical models for the electric con- 
tribution to it. Here, we show that in the interaction of graphene described by the Dirac 
model with metallic plate the results obtained are not sensitive to the approach used (ei- 
ther Drude or plasma). This is, however, not the case when graphene is described by the 
hydrodynamic model (see Sec. V). 

A. Free energy as a function of temperature 

Numerical computations of the Casimir free energy per unit area between graphene de- 
scribed by the Dirac model and Au plate were performed by using Eqs. (pQ), (jl])-© 
and (Q with fii = 1. The dielectric permittivity of Au along the imaginary frequency axis 
was described either by the generalized Drude-like model with temperature-dependent re- 
laxation parameter—^ or by the generalized plasma-like model.— iMt£i These models use the 
six-oscillator approximation for the optical data extrapolated to zero frequency by means of 
simple Drude and plasma models, respectively, with the plasma frequency u p = 9.0 eV and 
the relaxation parameter at room temperature 7 = 0.035 eV. At lower T the lower values of 
7 according to the standard theory of electron-phonon interaction have been used. 57 

The computational results using the Drude- and plasma-model approaches are found to be 
indistinguishable. Thus, at T = 77 K the relative difference between the Casimir free energies 
computed using both approaches achieves the maximum value of 0.02% at a = 50 nm, does 
not depend on A in the limits of our computational accuracy, and decreases with the increase 
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of separation. At T = 300 K this difference achieves the maximum values of 0.06% at A = 
and 0.07% at A = 0.1 eV, and again decreases with increasing a. Figure [7] presents the 
Casimir free energy per unit function of temperature (a) for a = 100 nm and (b) 

for a = 1 /mi. It can be seen that Fig. [7J demonstrates the same characteristic features, 
as Figs. [T] and H] plotted for dielectric plates, but the magnitudes of the free energy for 
the case of metallic plate are larger. The most important novel qualitative effect, which is 
found for both dielectrics and metals, is that for each nonzero A the free energy is nearly 
unchanged with the increase of T within some temperature interval. In Fig. [8] we present 
the computational results for the absolute thermal correction to the Casimir energy at zero 
temperature, defined in Eq. (1101) . as a function of temperature. The lines 1, 2, 3, and 4 
correspond to the values of the mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, 
respectively. This figure is analogous to Fig. [2] plotted for Si. It demonstrates that for a 
metallic plate the thermal correction in the graphene-plate geometry behaves qualitatively 
in the same way as for a dielectric plate, but with slightly larger magnitudes of the thermal 
correction. 

The computational results for the relative thermal correction, defined in Eq. ( TTTi) are 
presented in Fig. [9] as a function of temperature. The lines 1, 2, 3, and 4 again correspond 
to the same respective A, as in Fig. [HJ Figure [9] is analogous to Fig. [3] plotted for a dielectric 
plate (silicon). For a metallic plate at a = 100 nm the relative thermal correction at room 
temperature appears only slightly larger than for a dielectric plate. At a = 1 fim at room 
temperature the relative thermal correction for Au is smaller than for Si. This is explained 
by different values of the Casimir energy at zero temperature. 

B. Free energy as a function of separation 

Keeping in mind that in most experiments on the Casimir force the temperature is pre- 
served constant and measurements are performed at different separation distances, here 
we present the computational results for a free energy of graphene-metal interaction as a 
function of separation. In Fig. ITOTa) the free energy of graphene interacting with Au plate 
normalized on the Casimir energy between ideal metal planes ( ff2l) is shown. The three lines 
from bottom to top correspond to the values of mass gap parameter A = 0.1, 0.05, and 
< 0.01 eV, respectively. The obtained values of the free energy are larger than for graphene 
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interacting with Si plate [compare with Fig. Efa)]. At T = 300 K the dependence of the com- 
putational results on the mass gap parameter becomes not so pronounced as in Fig. ITDT a). 
This can be seen in Fig. [TOT b) where the bottom and top lines correspond to A = 0.1 eV 
and < 0.01 eV, respectively. 

It is interesting to consider the Casimir interaction of graphene with a ferromagnetic 
metal. It was shown^^ that at room temperature the ferromagnetic properties of real met- 
als may influence the Casimir force only through the contribution of the zero-frequency term 
of the Lifshitz formula. Recently the gradient of the Casimir force between a nonmagnetic 
Au sphere and a magnetic metal (Ni) plate has been measured.— We have computed the 
Casimir free energy per unit area between a graphene described by the Dirac model and 
Ni plate using the same formalism, as for an Au plate. The dielectric permittivity of Ni 
along the imaginary frequency axis was found from the tabulated optical data^i extrapolated 
to zero frequency either by the Drude or by the plasma model with the plasma frequency 
Up = 4.89 eV and the relaxation parameter at room temperature 7 = 0.0436 eV.— ^ The 
value of /i(0) = 110 for the static magnetic permeability of Ni has been used. It was found 
that relative differences in the computational results for the free energy of graphene-Ni inter- 
action, when Ni is described using the Drude- and plasma-model approaches, are as small as 
computed above for the interaction of graphene with an Au plate. The influence of magnetic 
properties on the free energy was also shown to be negligibly small. The relative difference 
between the free energies of graphene-Ni and graphene- Au interactions for graphene with 
A = 0.1 eV computed at T = 300 K is equal to 6% at a = 100 nm and decreases to 1% at 
a = 1 fim. It is less for smaller values of the mass gap parameter. Note that even these small 
differences are not due to magnetic properties of Ni but due to different plasma frequency 
and optical properties of Ni as compared to Au. 

C. Asymptotic behavior at high temperature 

Now we derive the analytic expression for the Casimir free energy of graphene described 
by the Dirac model interacting with metallic plate at r ^> 1. The contribution of the TM 
reflection coefficient for graphene interacting with dielectric plate was obtained in Eq. (l22l) . 
Taking into account that for metallic materials ro defined in Eq. (1201 is equal to unity, one 
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obtains from Eq. ( 122]) 



^tm(«, T) 



k B T 



ydy 



ln(l - e- y ) + - 
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n 00 (o) i - e- 



(24) 



Calculating the integral with respect to ?/ and using Eq. (fTTj) . we arrive at 
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where C(z) is the Riemann zeta function. This result in the special case A = was obtained 
in Ref.— . Note that for a dielectric plate considered in Sec. IIIC the contribution of the TM 
mode was in fact equal to the total free energy because r^(0,y) = for dielectrics. For 
metals this is in general not so (see below). 

The TE reflection coefficient for graphene at zero Matsubara frequency is obtained from 
Eq. © 

Mfn ^ n tr (0,y)-rioo(0,y) 



n tr (0,y)-Il 00 (0,y) + 2y 
From Eq. ([6]) taken at I = 0, Co = and Eq. f lT3|) it is easily seen that 

x{l — x) 
1 



(26) 



Tl tr (0, y) - n 00 (0, y) = 8av F y / dx 



x 1 



1 + e- 



1 + cosh ^9 



(27) 



where the quantity 9 is defined in Eq. (1141) . After identical transformations with account of 
Eq. (fTBl the result is 

1 , x(l-x) , n9 
■' - — tanh — 

9 T 
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tanh . 
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(28) 



This quantity is negligibly small as compared to unity because the main contribution to the 
Lifshitz formula (CQ) is given by y ~ 1 and for A — > one has 

Aiiav 2 F y 2 
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(29) 
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Thus, we can neglect by the difference of polarization operators in the denominator of 
Eq. (T26D and get 

r«(0,,)«-^tanh^. (30) 

Using the Lifshitz formula ([1]) and Eq. ( 130|) for the contribution of the TE reflection coefficient 
to the Casimir free energy of graphene-metal interaction at high temperature, we arrive at 

1 rji poo 

Jr E (a,T)^-^ ydyki 1 - rg(0,j/)rg(0,2/) e -«' 
Ibna' Jo L 



k B T av 



,~2 
F 



2A-KO? A 



A POO 

tanh _ / y *dyr®{0,y)e-*. (31) 
r Jo 



Here we have used that \r^(0, y)\ <C 1 at y ~ 1. 

Now we are in a position to consider metallic plates made of nonmagnetic and magnetic 
metals described within both the Drude and the plasma model approaches and in all cases 
find the high-temperature behavior of the total Casimir free energy. We begin with a non- 
magnetic metal described by the Drude model approach. In this case from Eq. ([9]) one 
obtains that r^(0,y) = and in accordance with Eq. f[3~Tj) the TE contribution to the free 
energy J-te^T") vanishes. Thus, for the plate made of a nonmagnetic Drude metal the 
total free energy of graphene-metal interaction at high temperature is given by Eq. (1251) . 

Next we consider a nonmagnetic metal described by the plasma-model approach. In this 
case from Eq. fl9]) we have 



r«(0, y) = g - V^±g! 2*V), (32) 



5y + y/1 + 6 2 y 2 
where the parameter 5 is defined as 

, s «« = « 4 <t (33) 

oo p zaujp 2a 

and 5 P = c/u p is the effective penetration depth of electromagnetic oscillations into the 
metal. Substituting Eq. (T32|) in Eq. ( 131]) and integrating with respect to y, one obtains 

Jr E (a,T) « -M-^Sitanh— (1 -G6 + 246 2 ). (34) 
12ira z A t 

The total asymptotic expression for the free energy at high temperature is given by the 
sum of ([25]) and (!34j) . Note that the contribution of the TE mode (!34j) is a negligibly small 
correction because av\ ~ 10 -7 . 
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We are coming now to the consideration of magnetic metals described by the Drude 
model. From Eq. (Q it follows 

rffi«U) = gg. (35) 
The substitution of this reflection coefficient in Eq. (T3TT) results in 

k B T avl vrA m(0) - 1 f°°2 -vj 

^ TE(a ' T) K «x tanh -r Moyn /„ y 6 dy 

- fc * T ^tanh^^. (36) 



127ra 2 A r /i(0) + l' 

This term is again negligibly small, as compared to J-tm(o, T), so that the total Casimir 
free energy at high temperature is well described by Eq. ( 125] ). 

Finally we consider the asymptotic expression for the free energy of graphene interacting 
with a magnetic metal described by the plasma model. In this case from Eq. (Q we get 



gto.,)- '^-^ 

1 - 25^/JI(0)y + 26 2 fi(0)y 2 

This is similar to Eq. ( 132]) with the replacement of 5 for 5 \J /x(0). Thus, instead of Eq. (|34|) . 
one obtains 

it i' r r\ ~ fcBT a ^ 4- v. n ^ 

.Ate (a, T) « -— — - tanh 

12na z A t 



x 



1 - 65v/^(0) + 24<J 2 /x(0)l • (38) 



The total asymptotic expression for the free energy at high temperature is given by the 
sum of Eq. (125]) and negligibly small addition (138]) depending on the properties of magnetic 
metal. The obtained analytic expressions were found in good agreement with the results of 
numerical computations within appropriate temperature (separation) intervals. 



V. COMPARISON BETWEEN HYDRODYNAMIC AND DIRAC MODELS OF 
GRAPHENE 

On this section we compare the computational results for the free energy of graphene- 
plate interaction obtained using two different models of graphene discussed in Sees. I and II. 
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We find separation regions where the predictions of both models are distinct and similar and 
compare respective asymptotic expressions for the free energy at high temperature (large 
separations) . 

A. Comparison between computational results for graphene described by two dif- 
ferent models 

Keeping in mind the possibility to compare theoretical predictions for the Casimir force 
with the experimental data, we calculate the free energy of graphene-plate interaction as a 
function of separation for both dielectric and metallic plates. Computations were performed 
using the Lifshitz formula (JT]) where the reflection coefficients (ITT)) for graphene in the 
framework of the hydrodynamic model were used. The computational results for T '/ 'Eq as 
a function of separation are presented in Fig. [TT] by the dashed lines (a) at T = 77 K and 
(b) at T = 300 K. In the same figure the respective results for T jE c computed using the 
Dirac model of graphene are reproduced from Fig. |5^a,b) by the solid lines. The solid lines 
in Fig. ITTT a) from bottom to top correspond to A = 0.1, 0.05, and < 0.01 eV, respectively. 
Note that in the scale used in Fig. ITTT b) the two lines of Fig. E^b) overlap. They are shown 
as a single solid line in Fig. ITTT b) . 

As is seen in Fig. ITTTa). at T = 77 K the hydrodynamic model of graphene predicts much 
larger magnitudes of the Casimir free energy than the Dirac model. Thus, at T = 77 K the 
predictions of the hydrodynamic model for \T\ at a = 0.5 /zm is by factors of 36.0 and 6.7 
larger than predictions of the Dirac model with A = 0.1 eV and A < 0.01 eV, respectively. 
At a = 1.5 /im the respective factors are 97.9 and 4.2. At T = 300 K [see Fig. ITTT b)] the 
predictions of the hydrodynamic model are larger than the predictions of the Dirac model 
by the factors of 2.5 and 1.3 at a = 0.5 /im and a = 1.5/im, respectively. As is seen in 
Fig. ITTT b). at T = 300 K, a > 4//m the asymptotic regime of large r is already achieved 
and the predictions of the hydrodynamic and Dirac models almost coincide. At T = 77 K 
[Fig. ITTT a)] the asymptotic regime of large r is achieved at much larger separations than 
those shown in the figure. 

Now we compare the predictions of the hydrodynamic and Dirac models of graphene 
interacting with a metallic plate. All computations were performed using the same formalism 
as above. We considered the plates made of a nonmagnetic metal Au and a magnetic 
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metal Ni. Each of these metals was described either using the Drude- or the plasma-model 
approach. 

The computational results for the normalized free energy J 7 / Ec are presented in Fig. [T5] 
at T = 300 K as a function of separation. In this figure, the solid line is reproduced from 
Fig. [TUT d) [note that in the scale of Fig. H2] the two solid lines in Fig. ITUT b) overlap] . This 
line is obtained using the Dirac model of graphene interacting with a metallic plate, be it Au 
or Ni (see Sec. IVB). The dashed and dotted lines 1, 2, 3, and 4 are obtained from compu- 
tations using the hydrodynamic model of graphene. The computational data shown by the 
dashed lines 1 and 3 were found for an Au plate described by the plasma- and Drude-model 
approaches, respectively. The data shown by the dotted lines 2 and 4 were computed for a 
Ni plate also described by the respective plasma- and Drude-model approaches. As can be 
seen in Fig. [121 the free energy of graphene-metal interaction with graphene described by 
the hydrodynamic model strongly depends on the metal used (magnetic or nonmagnetic) 
and on the chosen approach to the description of dielectric properties. Only for the non- 
magnetic metal (Au) described by the Drude model (line 3) the behavior of the free energy 
is qualitatively similar to the case of dielectric plate and nears the prediction of the Dirac 
model at large separations [compare with Fig. [TTT b)]. 

As an example, the Casimir free energy of graphene-metal interaction computed using the 
hydrodynamic model of graphene at a = 1 fim, T = 300 K is larger than the same quantity 
computed using the Dirac model by factors of 2.37 (for an Au plate described by the plasma 
model), 2.15 (for a Ni plate described by the plasma model), 1.93 (for an Au plate described 
by the Drude model), and 1.55 (for a Ni plate described by the Drude model). This allows 
comparison between different theoretical predictions and experimental data. It is interesting 
that for lines 1-3 the magnitude of J 7 predicted by the hydrodynamic model is always larger 
than for the predictions of the Dirac model (solid line). As to the line 4 (Ni described by the 
Drude model), the prediction for \JF\ from the hydrodynamic model becomes less than from 
the Dirac model at a w 1.55 /im and remains so at larger separations. Thus, at a = 5 fim 
the ratio between the predictions of hydrodynamic and Dirac models is equal to 0.42. 
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B. Asymptotic behavior at high temperature 



As was mentioned in Sec. IIIC, at high temperature (large separations) the zero-frequency 
term of the Lifshitz formula (prj alone determines the Casimir free energy. The reflection 
coefficients on the graphene, described by the hydrodynamic model, at zero frequency follow 
from Eq. (J7J) 

rS(0,y) = l, rg(0,y) = -^-. (39) 

K + y 

Taking into account that for a hydrodynamic model at T = 300 K the high-temperature 
regime starts at a > 5 /im, we find from Eq. (jSJ) that at these separations K > 6.75. 
Substituting Eq. (139]) for the reflection coefficient r^M info the zero-frequency term of the 
Lifshitz formula (JI]), for the TM contribution to the free energy at high temperature one 
obtains 

t, rp poo 

jT TM (a,T) = -^ ydy\n 1 - rg,(0, y)e^ . (40) 
lD7ra J i J 

For a dielectric plate, using Eq. (120]) and integrating in Eq. (f40|) . we get 

J- TM (a,T) = -_^_Li 3 (r ). (41) 

As noted in Sec. IIIC, for a dielectric plate the contribution of the TE mode to the free 
energy vanishes. Thus, Eq. ( 14 ip provides the complete expression for the free energy of 
graphene-dielectric interaction at high temperature. Equation (141 p coincides with the first 
term in Eq. ( 1231) obtained for graphene described by the Dirac model [remind that for the 
Dirac model of graphene the asymptotic expression ( 1231) becomes applicable at much smaller 
separations at the same room temperature; see Fig. [TTT b)]. 
For a metallic plate, replacing ro with unity, we arrive at 

J- TM (a,T) = -^C(3). (42) 

This result coincides with the first term of the asymptotic expression ( 125]) obtained for 
graphene described by the Dirac model. Equation (142!) provides the total asymptotic ex- 
pression for the free energy only in the case when a nonmagnetic metal of the plate is 
described by the Drude model (see the dashed line 3 approaching the solid line in Fig. H~2l 
when the separation distance increases). 
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Now we consider the contribution of the TE mode to the free energy of graphene-metal 
Casimir interaction. From Eq. f )39|) we have 

1 
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(43) 
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where /3 = 1/K takes the maximum value /3 max ~ 0.15 and decreases with further increase of 
separation. Then, for the contribution of the TE mode to the free energy at high temperature 
one obtains 
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For a nonmagnetic metal described by the plasma model the reflection coefficient r^(0, y) 
is given by Eq. ([32]). Substituting Eq. ([32]) in Eq. we find 
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Taking into account that at a > 5 /im it holds /3 3> 45, a more simple expression is also valid 
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By combining Eq. (I42p and Eq. (1461) . the total free energy for the interaction of graphene 
with a nonmagnetic metal described by the plasma model is obtained 



F(a,T) 
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C(3) (1 - /3 + 1.184/3 2 ) 



(47) 



Note that the main contribution to this free energy is twice that in Eq. ( 12"5"]) related to 
graphene described by the Dirac model. This explains different behaviors of the dashed line 
labeled 1 and the solid line at the largest separations in Fig. [T2j 

We are coming now to a magnetic metal described by the Drude model. In this case the 
reflection coefficient of the plate at zero frequency is given by Eq. ( 135]) . The substitution of 
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Eq. (J35J) in Eq. 



leads to 
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After the integration in Eq. (l4*8|) the result is 
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3/3 2 [Li 3 (-r M ) + Li 4 (-r M )]}. 



(50) 

By combining Eq. (j42p and Eq. (150]) . we obtain the following total free energy for graphene 
interacting with a magnetic metal described by the Drude model: 
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As an example, for Ni = 0.982 and Eq. ( 1511) takes the form 
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In is seen that the main contribution to this expression differs from the main term in Eq. ( 125]) 
obtained for the Dirac model of graphene. This is reflected also in Fig. [12] (compare the 
dotted line labeled 4 and the solid line). 

For a magnetic metal described by the plasma model we use the reflection coefficient (137]) 
and substitute it in Eq. ( 144"]) . All calculations are similar to the case of nonmagnetic metal, 



but the quantity 5 is replaced for 5^ /j,(0). Thus, instead of Eq. ( |45l) . we obtain 
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In this equation, however, it is impermissible to neglect by the quantity 5\J /x(0) as compared 
to (3. Substituting the numerical values of constants to Eq. ( 153]) . we find 
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By combining Eq. ( )54|) and Eq. (142]) . we arrive at the total free energy of graphene interacting 
with a magnetic metal described by the plasma model 



The main contribution to Eq. (155]) is by a multiple two larger than the main contribution to 
Eq. ( 125]) obtained for the Dirac model of graphene (compare the dotted line labeled 2 and 
the solid line in Fig. [T2]) . 

The obtained analytic asymptotic expressions for the free energy of graphene described 
using the hydrodynamic model in graphene-metallic plate geometry is in good agreement 
with the results of numerical computations. Thus, for an Au plate described by the Drude- 
and plasma-model approaches at a = 5 /im, T = 300 K, the results of analytic and numerical 
calculations differ by 0.35% and 1.4%, respectively. For Ni the same relative differences are 
equal to 9.3% and 3%. Note that relatively large deviation obtained for Ni plate described 
by the Drude-model approach is explained by the fact that in this case at a = 5 /im the 
linear asymptotic regime is not yet achieved (see the dotted line labeled 4 in Fig. [T2]) . 

To conclude, we emphasize that all results for the Casimir free energy obtained in this and 
previous sections are simply convertable to the Casimir force F(a,T) in the experimentally 
relevant configuration of a sphere above a plate used in most of experiments on measuring 
the Casimir force. This can be done by means of the PFA which states that 



where R is the radius of the sphere. In our case, where J-"(a, T) is the Casimir free energy 
between a graphene sheet and a material plate, the force F(a, T) defined in Eq. fl56|) can 
be considered as the Casimir force acting between a graphene sheet and a material sphere 
of radius R. As was mentioned in Sec. I, the error introduced by the use of the PFA was 
recently proved^ - — to be smaller than a/R (i.e. or order of 0.1% for the typical values of 
parameters). Keeping in mind that many effects considered above far exceed 100%, the use 
of the PFA in the comparison between experiment and theory is fully justified. 




(55) 



F(a,T) = 2TrRT(a,T), 



(56) 
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VI. CONCLUSIONS AND DISCUSSION 



In the foregoing, we have investigated the Casimir free energy and the thermal correction 
to the Casimir energy at zero temperature for a suspended graphene sheet interacting with a 
material plate, either dielectric or metallic. In so doing graphene was described by the fully 
relativistic Dirac model with temperature-dependent polarization tensor. The dielectric 
properties of the plate were described by the frequency- dependent dielectric permittivity 
taking into account the interband transitions of core electrons. For a metallic plate both the 
Drude- and plasma-model approaches suggested in the literature have been used. 

The main novel result obtained for both dielectric and metallic plates is that for graphene 
with any nonzero mass gap parameter A there exists temperature interval where the Casimir 
free energy remains nearly constant. This happens under the condition ksT <C A, which 
should be satisfied with a large safety margin. If this condition is satisfied, the thermal 
correction to the Casimir energy at zero temperature remains negligibly small. We have also 
demonstrated that under the condition A < /c B T the thermal correction becomes relatively 
large. This makes possible large thermal corrections for a graphene sheet interacting with 
material plate at rather low temperature (short separations). 

With respect to the interaction with a metallic plate, it was shown that for graphene 
described by the Dirac model the computational results for the free energy are nearly in- 
dependent on whether the Drude- or plasma-model approach to the dielectric permittivity 
of metal is used. To a large extent the free energy of graphene interacting with metallic 
plate is also independent on whether metal is nonmagnetic or magnetic if graphene is de- 
scribed by the Dirac model. In all cases considered (dielectric or metallic plate, nonmagnetic 
or magnetic, described by the Drude- or plasma-model approach) the analytic asymptotic 
expressions for the Casimir free energy at high temperature (large separations) have been 
obtained and compared with the results of numerical computations. 

The Casimir free energies obtained using the Dirac model of graphene were compared 
with those calculated using the hydrodynamic model. It was shown that at moderate tem- 
peratures (separations) the magnitudes of the free energy computed using the hydrodynamic 
model of graphene differ significantly from that computed using the Dirac model. This can 
be used for the experimental test of these models. At large separations (high temeperature) 
the theoretical predictions from both models of graphene nearly coincide for a dielectric plate 
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and for a nonmagnetic metallic plate described by the Drude model. For a nonmagnetic 
metallic plate described by the plasma model and for a magnetic plate described by any 
model the hydrodynamic and Dirac descriptions of graphene lead to quite different results 
for the free energy at large separations (high temperature). This fact can be also used for 
the experimental test of different models. In this respect the investigation of the interaction 
between graphene and metamaterials^ 4 - is also of large interest. We have also found analytic 
asymptotic expressions for the free energy at high temperature (large separations) when 
graphene described by the hydrodynamic model interacts with a dielectric plate or with a 
plate made of a nonmagnetic or magnetic metal. In the last two cases both the Drude- 
and plasma-model approaches have been used for a description of the dielectric properties 
of metal. The calculation results obtained from the asymptotic expressions were found in a 
very good agreement with the results of numerical computations. 
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FIG. 1: (Color online) The Casimir free energy per unit area as a function of temperature for the 
interaction of graphene described by the Dirac model with Si plate (a) at a = 100 nm and (b) at 
a = 1 //m. The lines from top to bottom correspond to the mass gap parameter A = 0.1, 0.05, 
0.01, and < 0.001 eV, respectively. 
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FIG. 2: (Color online) The thermal correction to the Casimir energy per unit area at zero temper- 
ature as a function of temperature for the interaction of graphene described by the Dirac model 
with Si plate (a) at a = 100 nm and (b) at a = 1 /im. The lines labeled 1, 2, 3, and 4 correspond 
to the mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, respectively. 
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FIG. 3: (Color online) The relative thermal correction to the Casimir energy at zero temperature 
as a function of temperature for the interaction of graphene described by the Dirac model with Si 
plate (a) at a = 100 nm and (b) at a = 1 ^m. The lines labeled 1, 2, 3, and 4 correspond to the 
mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, respectively. 
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FIG. 4: (Color online) The Casimir free energy per unit area as a function of temperature at 
a = 100 nm for the interaction of graphene described by the Dirac model with sapphire plate. 
The lines from top to bottom correspond to the mass gap parameter (a) A = 0.1, 0.05, 0.01, and 
< 0.001 eV and (b) A = 0.1, 0.09, 0.08, 0.07, 0.06, and 0.05 eV, respectively. 
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FIG. 5: (Color online) The Casimir free energy as a function of separation for the interaction 
of graphene described by the Dirac model with Si plate (a) at T = 77 K and (b) at T = 300 K 
normalized on the Casimir energy between two parallel ideal metal planes. The lines from bottom 
to top correspond to the mass gap parameter (a) A = 0.1, 0.05, < 0.01 eV and (b) A = 0.1, 
< 0.01 eV, respectively. 
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FIG. 6: (Color online) The Casimir free energy as a function of separation for the interaction of 
graphene described by the Dirac model with sapphire plate at T = 77 K normalized on the Casimir 
energy between two parallel ideal metal planes. The lines from bottom to top correspond to the 
mass gap parameter (a) A = 0.1, 0.05, < 0.01 eV and (b) A = 0.1, 0.09, 0.08, 0.07, 0.06, 0.05 eV, 
respectively. 
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FIG. 7: (Color online) The Casimir free energy per unit area as a function of temperature for the 
interaction of graphene described by the Dirac model with Au plate (a) at a = 100 nm and (b) at 
a = 1 fim. The lines from top to bottom correspond to the mass gap parameter A = 0.1, 0.05, 
0.01, and < 0.001 eV, respectively. 
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FIG. 8: (Color online) The thermal correction to the Casimir energy per unit area at zero temper- 
ature as a function of temperature for the interaction of graphene described by the Dirac model 
with Au plate (a) at a = 100 nm and (b) at a = 1 ^m. The lines labeled 1, 2, 3, and 4 correspond 
to the mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, respectively. 
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FIG. 9: (Color online) The relative thermal correction to the Casimir energy at zero temperature 
as a function of temperature for the interaction of graphene described by the Dirac model with Au 
plate (a) at a = 100 nm and (b) at a = 1 ^m. The lines labeled 1, 2, 3, and 4 correspond to the 
mass gap parameter A = 0.1, 0.05, 0.01, and < 0.001 eV, respectively. 
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FIG. 10: (Color online) The Casimir free energy as a function of separation for the interaction 
of graphene described by the Dirac model with Au plate (a) at T = 77 K and (b) at T = 300 K 
normalized on the Casimir energy between two parallel ideal metal planes. The lines from bottom 
to top correspond to the mass gap parameter (a) A = 0.1, 0.05, < 0.01 eV and (b) A = 0.1, 
< 0.01 eV, respectively. 
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FIG. 11: (Color online) The Casimir free energy as a function of separation for the interaction of 
graphene described by the hydrodynamic model (dashed lines) and by the Dirac model (solid lines) 
with Si plate (a) at T = 77 K and (b) at T = 300 K normalized on the Casimir energy between 
two parallel ideal metal planes. The solid lines correspond (a) from bottom to top to the mass gap 
parameter A = 0.1, 0.05, < 0.01 eV, respectively, and (b) to A = 0.1 eV. 
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FIG. 12: (Color online) The Casimir free energy as a function of separation for the interaction 
of graphene described by the hydrodynamic model (dashed and dotted lines) and by the Dirac 
model (solid line) with Au (dashed lines labeled 1 and 3) and Ni (dotted lines labeled 2 and 4) at 
T = 300 K normalized on the Casimir energy between two parallel ideal metal planes. For lines 
labeled 1 and 2 metal is described by the plasma-model approach, and for lines labeled 3 and 4 by 
the Drude-model approach. 
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